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Abstract 

We study the support (i.e. the set of visited sites) of a t step random 
walk on a two-dimensional square lattice in the large t limit. A broad 
class of global properties M(t) of the support is considered, including, 
e.g., the number S(t) of its sites; the length of its boundary; the number 
of islands of unvisited sites that it encloses; the number of such islands 
of given shape, size, and orientation; and the number of occurrences 
in space of specific local patterns of visited and unvisited sites. On 
a finite lattice we determine the scaling functions that describe the 
averages M(t) on appropriate lattice size dependent time scales. On 
an infinite lattice we first observe that the M{t) all increase with t 
as ~ tj log* 1 1, where k is an M dependent positive integer. We then 
consider the class of random processes constituted by the fluctuations 
around average AM(t). We show that to leading order as t gets large 
these fluctuations are all proportional to a single universal random 
process rj(t), normalized to n 2 (t) = 1. For t —> oo the probability 
law of rj(t) tends to that of Varadhan's renormalized local time of self- 
intersections. An implication is that in the long time limit all AM(t) 
are proportional to AS(t). 
PACS 05.40+j 
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1 Introduction and summary 



The reason for the multiple connections between the random walk and many 
questions of current and of permanent interest in science is evidently mathe- 
matical: The random walk models the action of the Laplace operator, which 
is common to all those problems. Indeed, important early results on the 
random walk are due to mathematicians. Among them, the famous Polya 
theorem asserts that in spatial dimensions d < 2 a random walk is certain 
to return to its initial position, whereas in d > 2 it will escape to infinity. 
The random walk in d = 2 is therefore at its critical dimension for return to 
the origin. Many features associated with this fact make the random walk 
in two dimensions particularly interesting. For example, the probability dis- 
tribution of the time interval To between two successive visits of the walk to 
its initial position decays for large intervals as ~ l/rolog 2 ro, so that very 
long excursions occur away from the point of departure. Excellent recent 
monographs on random walk theory have been written by Weiss [|| and by 
Hughes g. 

In this work we consider the simple random walk on a two-dimensional 
square lattice: The walk starts at time t = at the origin x = and steps at 
t = 1, 2, 3, . . . with equal probability to one of the four neighboring lattice 
sites. Our investigation focuses on the statistical properties of the support 
of the walk at time t, i.e. of the set of sites that have been visited during 
the first t steps. 

At any given time t the support can be visualized as a set of black sites 
(the visited ones) in a lattice of otherwise white (unvisited) sites, as shown in 
Fig. 1. The set of unvisited sites is divided into components that are disjoint 
(i.e. not connected via any nearest neighbor link) and that we call islands. 
In the course of time existing islands will be reduced in size and single-site 
islands will eventually be destroyed; the number of islands increases each 
time that a step of the walk cuts an existing island, or the outer region 
surrounding the support, into disjoint components. 

Questions about the statistical properties of the islands are natural 
and our interest in them was raised by a simulation study by Coutinho, 
Coutinho-Filho, Gomes, and Nemirovsky y|. These authors investigated 
the evolution of the number of islands I(t) and several of their properties on 
finite lattices of up to 1200 2 sites. In a short report § we showed that an 
analytic calculation is possible, both on finite and infinite lattices, for some 
(not all) of the quantities considered by Coutinho et al. 0j, and that good 
agreement between theory and computer simulation is obtained. In this 
work we present a full account of most of the results announced in Ref. ||] 
and consider a great many related questions. We are, in particular, led to 
consider the infinite lattice again. 

It appears that the number of islands I(t) is but one member of a much 
wider class of observables, generically to be denoted as M(t), with closely 
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related properties. This class includes the total number S(t) of sites in the 
support as well as the total length E(t) of the boundary of the support (on 
an infinite lattice the boundary length is the sum of the external perimeter 
and the perimeters of the islands enclosed). 



The main variable that characterizes the support is its total number of 
sites S(t), sometimes called its range, which was first studied by Dvoretzky 
and Erdds || . In the limit of large times t it has the average value || 0, ||, ||] 



S(t) 



irt 
log 8t 



and the root mean square deviation [10, 11 



1/2 



AS 2 (t) ' ~ A 



7TU 



log fc VU 



(1.1) 



(1.2) 



where A = oo.B/B... and where we write AS" = S — S. The typical support 
is known to be far from spherical, and its principal moments of inertia and 
asphericity have been studied [12, 0|. Expressions exist [14] also for its 



span, that is, the smallest rectangular box that it fits in. 



The question of calculating the average S(t) becomes, in general dimen- 
sion and in an appropriately taken continuum limit, the celebrated Wiener 
sausage problem: What is the volume swept out in time t by a (i-dimensional 
Brownian sphere of finite radius? The Wiener sausage appears naturally in 
certain applications of the random walk such as, for example, the study 
by Kac and Luttinger yl| of Bose-Einstein condensation in the presence of 
impurities. Historically this continuum problem precedes its lattice coun- 
terpart. It was considered as early as 1933 by Leontovitsh and Kolmogorov 
| |16| | and is still today an active subject of investigation in probability theory 
( |17| ; sec for a recent overwiew) and in mathematical physics JT^] . 

The islands in the support of a lattice random walk have their counter- 
part in the connected components into which a two-dimensional Brownian 
motion path divides the plane. These have been studied by Mountford ppj ] 
and Le Gall 17], and, most recently, by Werner [21], who determines, among 
other things, how many there are larger than a given size e, in the limit of 
e — > 0. The outer boundary of the Brownian motion path has been consid- 
ered very recently by Lawler [23]. 



We characterize as follows the class of observables M(t) to be studied 
below. For each lattice site x we introduce an occupation number m(x,t), 
equal to if site x is visited by the walker before or at time t, and equal to 
1 if it is not. Let now A\ and A2 be disjoint finite sets of lattice vectors. 
Then the product 

Y[ m(x + ai,t) ]J [1 - m(x + a 2 , t)] (1.3) 

aieAi a 2 GA 2 
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codes for a specific spatial pattern a = (A\,A2) of white (unvisited) and 
black (visited) sites. When x runs through the lattice, the product ( |1.3| ) 
equals 1 when the pattern a is encountered and otherwise; hence this 
product summed on all x represents the total number N a (t) of occurrences 
in space of the pattern a. The observables that this work deals with are 
these pattern numbers N a (t) and their linear combinations M{t). It is easy 
to see that suitably chosen M(t) may represent, e.g., the total number S(t) 
of sites in the support; the total boundary length E(t) of the support; the 
total number I(t) of islands of unvisited sites enclosed by the support; or 
the total number Ip(t) of islands of a given type (3 (where type stands for 
shape, size, and orientation). Our main results concerning the observables 
M{t) are of two kinds. 

(i) Average behavior on a finite lattice. On a finite lattice of N sites the 
averages M(t) approach their limiting values on the time scale t ~ iVlog 2 N, 
whereas their main variation occurs on the earlier time scale t ~ iVlogiV. 
We calculate the scaling functions that describe the time dependence on 
both time scales. Several of the observables M{t) of interest just mentioned 
are treated as examples. A comparison is made, where possible, with the 
simulations by Coutinho et al. § 

(ii) Average behavior and fluctuations on the infinite lattice. On an infi- 
nite lattice the fluctuating properties of the support manifest a universality 
that can be described as follows. Let M and M' be two linear combinations 
of pattern numbers. Then, with the same notation as before for deviations 
from average, we show by explicit calculation that asymptotically for t — > oo 

Kk+l t TT k ' + 1 t 

M(t) ~ m k rrn — , M'(t) ~ mi, rj—, — (1.4) 

W log fc+1 8t W h log k+1 8t V ' 

^ll+ll'+e u e 

wwrr?vn (L5) 



AM(t)AM'(t) ~ 4A e (|| + oo)(||' + oo 



Here k and k' are nonnegative integers that depend on M and M', respec- 
tively, and are called their order; and m' k , are known proportionality 
constants; and the same number A appears that was first encountered in 
the study of the variance of S(t). The above equations imply that the nor- 
malized deviations from average 

log8* AM jt) 

m{t) = (k+i)A m ( } 



have a correlation matrix r]M{t) r ]M l {t) whose elements all equal unity. This 
can be true only if all t]m are equal to a single random variable to be called 
rj{t). As a consequence the random variables AM(t) are, to leading order 
as t — > oo, all proportional to n(t). Explicitly, 



7T k+1 t 



AM(t) ~(k + l)m k k+2Q A V (U) (1.7) 
log' c+z 8t 
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where r](t) is universal (independent of M). Hence what seemed to be a 
large number of independent fluctuating degrees of freedom of the support 
is hereby reduced, to leading order as t — > oo, to only a single degree of 
freedom. 

We now briefly summarize the contents of the successive sections. Since 
all quantities of interest are, in the end, expressed in terms of the random 
walk Green function, we collect in Sec. ^| the basic knowledge required about 
this function. In Sec.^ we express the total boundary length E(t) and the 
total number of islands I(t) in terms of the occupation numbers and discuss 



the wider class of observables M(t) of which they are examples. In Sees. 4.1 



and [4.2| we calculate the time dependent averages M(t) of such observables. 



In Sees. 4.3-4.6 the general expression for the result, Eq. (4.16), is analyzed 



for large times and explicitly worked out for several examples both on the 



infinite and the finite lattice. In Sees. 5.1 and 5.2 we show, for the infinite 



lattice, how to calculate the correlation between two observables M(t) and 
M'(t). In Sec. |5.3| we arrive at the final results concerning the random vari- 
able r](t). In the discussion in Sec.^ we present various comments, compare 
where possible our lattice results to their continuum analogs, and speculate 
about further connections between several quantities. 

2 The random walk Green function 
2.1 Definition 

A random walker starts at time t = at the origin x = of a square lattice 
and steps at each instant of time t = 1, 2, 3, . . . with probability \ to one 
of its four nearest-neighbor sites. The Green function G(x, t) denotes the 
probability that at time t the walker is at lattice site x, and 

oo 

G(x,*) = 5>*G(x,*) (2.1) 
t=o 

its generating function. In this section we collect in concise form those prop- 
erties of the generating function that will be needed later. An elementary 
calculation gives, for a finite periodic lattice of L x L = N sites, 

Gtvi^Er^i 1 t ( 2 - 2 ) 

Jy q 1 — 2 Z { cos Ql + cos Q2) 
where q = (qi,q%) = 2tt{k\, Kq) / L with the running through the values 



0, 1, 2, . . . , L — 1. In the limit of an infinite lattice expression (2.2) becomes 



i i"K r-K g— iq-x 
G(x,z) = -—2 d Ql dq 2 - — (2.3) 

(27TJ Z J-ir J-ir 1 — ^(COS qi + COS q2) 

At each point in this work it will be clear whether we are discussing the finite 
or the infinite lattice; in some cases we shall denote corresponding quantities 
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in the two geometries by the same symbol, as for example in Eqs. (2.2) and 
( |2.3D , and not explicitly indicate their N dependence on a finite lattice. 



2.2 Expansion near z = 1 

The long time behavior of the physical quantities of interest is determined 
by the behavior of G(x,z) in the complex plane near z = 1. Expressions 
(2.2) and (p. 3D both have the property that for z — ► 1 the function G(x, z) 



diverges. In order to study this divergence it is convenient to write 

G(x,z)=G(0,z)-g(x,z) (2.4) 

where on the RHS the term g(x, z) contains all the x dependence and re- 
mains finite for z — ► 1. We shall discuss G(0, z) and g(x, z) separately. 

2.2.1 The function G{0,z) 



Finite lattice. Expression (|2.2| ) has a simple pole as a function of z when- 
ever one of the denominators inside the sum on q vanishes. This leads to a 
sequence of poles on the real axis for z > 1, of which the first one is located 
exactly at z = 1. The interval between two successive poles is of 0(N^ 1 ) 
and contains a zero. Upon setting x = in Eq. (|2.2| ) and expanding each 
term for small 1 — z one gets 



G(0, z) = — L - + a(iV) - oi(JV)(l - *) + 0((1 - z) 2 ) 



(2.5) 



in which the coefficients are functions of N that in the limit TV" — ► oo behave 
as H 

a(JV) = ilogciV + C'(iV- 1 ) 

ai (N) = Cl N + 0(logN) (2.6) 
with c = 1.8456... and ci = 0.06187... 



The expansion in Eq. (2.5) represents well the behavior of G(0,z) near 
the pole at z = 1, but is certainly not valid on approach of the next pole. It 
can be used, however, to determine the location of the zero z = zq between 



the first two poles. Upon solving Eq. (2.5) in successive orders for zq one 
finds 



z = l + 



1 



1-^ + 



Na(N) L a 2 (N) 



(2.7) 



From Eqs. fl2,7|) and (2.6) we conclude that when N —> oo this zero is sepa- 
rated from the pole at z = 1 by a distance only of 0(1/N log N). All other 
zeros of G are separated from z = 1 by a distance of at least 0(N^ 1 ). Due 
to its exceptional proximity to z = 1 the zero zq plays a special role in the 
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long-time behavior of the random walk on finite lattices. This fact was first 
noted by Weiss, Havlin, and Bunde |Q| and has also been exploited |^] in 
the study of the covering time of a finite lattice by a random walk. 
Infinite lattice. When N — > oo the poles of G density to a branch cut and 
the expansion near z = 1 is [^] 

6(0, z)= X - log +0{(1- z) log(l - z)) (2.8) 

Later on in this work we shall also for brevity denote (5(0, z) as Gq(z). 



2.2.2 The function g(x,z) 

Expanding g(x, z) for finite iV around z = 1 gives 

<?(x, z) = g N (x.) + <^(x)(l -z) + 0((1 - zf) (2.9) 

where we now explicitly indicate that the expansion coefficients are N de- 
pendent. We shall set 

5 (x) = Jim g N (x). (2.10) 



Spitzer [27] shows how to calculate the g(x) for x close to the origin. Letting 
ei and e 2 denote the unit vectors we have the following values: 

0(0) = s(2ei) = 4- 8/vr 

g(e x ) = l g{2e l + e 2 ) = S/Ti-l (2.11) 

g(ei + e 2 ) = 4/vr 

When combining preceding results we see that for finite N the quantity 
G(x, z) has the expansion 

G(x, z) ~ — j- + [o(JV) - 5 ^(x)] + 0(1 - z) (2.12) 

whose second term behaves for large N as 

a(N) - g N (x) ~ 1 log ciV - 5 (x) + • • • (2.13) 

7T 

with the dots representing terms that vanish as ./V — > oo. 



2.3 Scaling limit 

In our study of the fluctuations in Sec.|| we shall also need G(x, z) in the 
scaling limit z — ► 1, x — > oo with x 2 (l — z) fixed. The behavior in this limit 

is! 

G(x, z) ~ -K (2x(l - z) 1 ' 2 ) (2.14) 
where Kq is the modified Bessel function of order zero. As is well-known, 



for z — * 1 the dominant contribution to the sum (2.1) comes from values 
of t that are of 0((1 — z)^ 1 ), and therefore the scaling limit corresponds to 
focusing on distances x of 0(y/i). 
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3 Islands and other observables 



3.1 Islands 

The representation of Fig. 1 divides the lattice into a black area and white 
ones. We introduce for each lattice site x the occupation number 

, . _ f 1 if x has not yet been visited at time t (3 1) 

|0 otherwise 

At this point it may seem more natural to work with the equivalent occu- 
pation numbers n(x, t) = 1 — m(x, t), but definition (|3.l| ) will soon prove to 
be more convenient. 

If starting at an arbitrary element, one follows the boundary such that 
the white sites are on the left and the black ones on the right, then one 
will return to the point of departure after having turned either through an 
angle 2ir (if that part of the boundary encloses an island in the support) 
or through an angle — 2ir (if it is the outer boundary of the support). The 
number of islands is therefore obtained from the number of turns in the 
boundary, by adding those of Fig. 2a with weight +j and those of Fig. 2b 
with weig ht -\. The dia grams of Fig. 2c correspond to two turns of n/2 
and therefore count with weight factor ^. This procedure counts the outer 
boundary with weight —1 and therefore 2 has to be added to obtain the 
final result. It is now easy to construct the expression for the total number 
of islands in terms of the occupation variables m(x,t). We have to shift a 
2x2 window across the lattice, check all 2 x 2 local site configurations and 
add all those that are of the types of Fig. 2 with the proper weights. Let x 
denote the lower lefthand site in the 2x2 window. Then, for example, the 
expression 

m(x, t) [1 — m(x + ei, t)] [1 — m(x + e2, t)] m(x + ei + e2, t) 

(3.2) 



(which is of type ( |1.3| )) equals 1 or according to whether the local 2x2 
configuration is or is not equal to the diagram of Fig. 2c. Writing down anal- 
ogous expressions for all other diagrams, summing these, and subsequently 
summing them on x yields I(t) as a linear combination of pattern numbers. 
After rearranging terms one finds 

I(t) = 2 + ^2 [f&( x ! t) ~ rn(~K, £)m(x + ei, t) — m(x, i)m(x + e2, t) 

X 

+ m(x, t)m(x + ei,i)m(x + e2, t)m(x + ei + e2, t) (3-3) 

This expression is at the basis of all calculations concerning islands on the 
infinite lattice. Finite lattice calculations require a further remark, which is 
made below Eq. (4.44). 
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3.2 Other observables 



We are now interested in other variables whose values can be obtained by 
the "window" method. To see the general form of these variables, let A 
be a finite subset of lattice vectors. The set {x + a| a G A}, obtained by 
translating A by a vector x, will be written as x + A. The variable 

m x+A (t) = J] m(x + a, t) (3.4) 

aeA 

is equal to unity if at time t all the sites of this set are white (unvisited), 
and is zero otherwise. The sum variable 

M A (t) = Y,m x+A {t) (3.5) 

counts the total number of wholly white sets in the lattice that can be 
obtained from A by a translation. 

In the remainder we shall also wish to take for A the empty set 0. In 
that case the right hand side of Eq. ( |3.4[ ) should be assigned the value unity 
and Eq. (|3.5| ) shows that M$(t) is the total number of lattice sites. 

The number of islands I(t) was initially found in Sec. |3. 1| as a linear 
combination of pattern numbers, which led to the final expression ( |3.3| ). The 
representation of an observable as a linear combination of pattern numbers 
is in general nonunique, but expression ( |3.3| ) is a unique member of the class 
of variables M (t) that are of the form 

M(t)=J2^AM A (t) (3.6) 
A 

with arbitrary coefficients fx A . Various other quantities of potential interest 
can be expressed this way. Some of these, with their coefficients fi A , have 
been listed in Table I. The best known example is the total number S(t) of 
sites in the support, 

S(t) = £[l-m(x,t)] (3.7) 

X 

which has \i A = ±1 for A = and A = {0}, respectively, and fi A = 
otherwise. Another example is the total boundary length E(t) between 
the visited and unvisited lattice sites (that is, the total number of pairs of 
neighboring sites of which one is white and one black). It can be expressed 
as 

E(t) = [m(x, t)(l - m(x + ei, t) + (1 - m(x, i))ra(x + e l5 1) 

X 

+ m(x, t)(l - m(x + e 2) t) + (1 — m(x, t))m(x + e 2 , t) 
= ^2 [4m(x, t) - 2m(x, i)m(x + e 1 ,t) - 2m(x, t)m(x + e 2 , i) (3.8) 
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We note that one has the relation 



£/M = (3.9) 
A 

for the three observables S, E, and /, but that 

E^ = (3-10) 

A^% 

only for E and I, but not for S. Property ( fj.Sj ) is required if the sum on x 
in Eq. ( |3.5| ) is to have a finite limit when the lattice size iV tends to infinity. 
The property ( 3.10| ) is easily traced back to the fact that for E and I the 



window selects patterns consisting of both white and black sites, whereas for 
S it selects only black sites. When Eq. ( |Q| ) holds, Eq. ( 3. 1C| ) is equivalent 
to /U0 = 0. 

4 Averages of observables 

4.1 Relation to first passage time probabilities 

We shall now calculate in a unified way averages and fluctuations of quan- 
tities M(t) of type ( |3.6| ). The approach of this subsection will serve as the 
basis for the developments to follow. Using Eqs. ( |3.4|) -( ^6[ ) we see that the 
averages M{t) are linear combinations of expressions of the type 



m x+A {t) = JJm(x + a,t) 

at A 



T=0 

= i-EEW x + a > T ) C 4 - 1 ) 

r=0 aeA 

in which the last two transformations make sense only when A ^ 0; / x +a( t ) 
is the probability that the walker's first visit to any site of the set x + A 
takes place at time r; and /x+^x + a, r) is the probability that it takes 
place at time r and that it concerns the specific site x + a. Upon averaging 
Eq. (|3.5|) , using Eq. ( |4.1| ), and passing to generating functions we find, for 
all nonempty A, 

^a(z) = e [ E W x + *>*) - x ] ( 4 - 2 ) 

For the empty set one derives directly that 

%(z) = N/(l - z) (4.3) 
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We now sum the M a{z) given by Eqs. fl4.2| ) and (4.3) on all A with coeffi- 
cients li A . Using Eq. ( j3li| ) and introducing for all A ^ 



F A (a,z) = J2h+A(^ + a,z) (4.4) 

X 

yields 

-M(z) = V-A E ^(a, z) (4.5) 

With this formula we have reduced the generating function of the average of 
interest, M(t), to the quantities Fa which are closely related to first passage 
times, but still unknown. We shall now proceed to determine the Fa- 



4.2 Solving the first passage time probabilities 

Standard random walk theory |^7|, ^] relates the first passage probabilities 
/ to the Green function G by 



G(x + a,z)= E / x+ a(x + a')G(a - a', z) (4.6) 

a'eA 

for all a G A With the aid of Eq. ( [4.4|) we find for ^4 the equation 

J2F A ( a ',z)G( a - a ',z) = -^— (4.7) 

for all a G A. This is a matrix equation for the Fa whose dimension is the 
number \A\ of sites in the set A. This equation possesses special properties 
which are best exhibited by converting it to the shorthand notation 

7aa'( 2 ) = 9(a-a.',z)/Go(z) 

F a = (1- z)G (z)F A (a,z) (4.8) 

Eq. (|4.7|) then becomes 



E( 1 -^a'(^))^"a' = l (4.9) 
a'ei 

for all a G A. If we denote by 7^ the matrix of elements 7 aa ' with a, a' G ^4, 
then in matrix notation 

(J- 7 (A) (^ = j (4.10) 

where J and j are the matrix and vector, respectively, of dimension \A\, 
whose elements all equal 1. As shown in Eq. (f4.5|), we only need the sum of 
the components of J 7 . Formal inversion gives 

aeA aeAa'eA 
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In a later stage we shall wish to take the limit z — ► 1. In view of Eq. 
and the known behavior of Gq(z) this implies that j( A '(z) — > 0, so that, 
except when |A| = 1, the matrix inverse (J — 7^(z)) _1 in Eq. ( [4.1 1| ) ceases 
to exist. We therefore now convert that equation to a form more suitable 
for taking that limit. In the appendix it is shown that 

S'-T^ (4 ' 12) 

where 

T^t^EE^wrV (4.i3) 

aeAa'eA 

Upon coming back to the original notation, but now with the abbreviation 

g A (z) = lA {z)G {z) (4.14) 
we find from Eq. ( |4. 12 ) the solution of Eq. ( [4.7|) in the form 



V F A (bl, z) = rrrr — z — r y-rr (4.15) 

£a (I - z)[G (z) - g A (z)} K ' 

When \A\ = 1, in which case we may take A = {0}, one deduces directly 



from ( |4 that ( 4.15 ) holds with g^(z) = 0. For A of diameter not too 
large, as is the case in many examples of interest, the quantity g A (z) is 
easily expressed explicitly in terms of the g(a — a', z). If after substitution 



of expression ( |4~T5l ) in Eq. Q we sum on A, use Eq. ( |3.9| ), and transform 



back to the time domain, we get 

where the integral runs counterclockwise around the origin. This result is 
still fully exact and applies to both finite and infinite lattices. 

4.3 Long-time behavior of averages. Infinite lattice 



The starting point for the analysis of this section is Eq. (4.16). On an infinite 
lattice the asymptotic behavior of M(t) as t — > oo is determined by the 
z — ► 1 behavior of the integrand. This behavior follows from expression 



(2.8) for Gq(z) and from Eqs. (|4.14p and ( [4.13 ) which together determine 



g A (z). We shall satisfy ourselves with retaining the leading z — > 1 behavior 
and corrections that are of relative order of negative powers of log(l— z). The 
terms neglected are of relative order 1 — z, apart from logarithmic factors. 
This means that in Eq. ( [4.16 ) we may replace Gq{z) with 7r _1 log(8/(l — z)) 



and the g A (z) with their values at z = 1, which for brevity we shall denote by 
g A . The coefficient g A appears in potential theory and is the two-dimensional 
lattice analog of the electrostatic capacity of the set A; its properties have 
been reviewed by Spitzer pTfl . 
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We expand the summand in Eq. ( 4.16 ) in inverse powers of log(8/(l— z)). 
Using the above results we so obtain, asymptotically for t — > oo, 



M(t) ~ m nJn+l,2{t) 
n=0 



(4.17) 



where the coefficients m n are determined by the observable M according to 

(4.18) 



m r . 



for n = 0, 1, 2, . . . and where 



Jni{t) = 



dz 



2m J (1 - z)^log n T ^ 5 



(4.19) 



The t — ► oo behavior of the integrals ( |4.19| ) has been studied, in particular, 
by Henyey and Seshadri |2£| for the case 1 = 2. A generalization of their 
result is 



Jnl{t) 



n,f—\ oo 
TT n t 



log™ 8t 



with coefficients 



m=0 



d" 



m + n — 1 
m 



log m 8t 



1 



cfe m r(x + i) 



x=0 



(4.20) 



(4.21) 



We have = !/(•£ — 1)! and shall also need explicitly below 



&12 
&13 
^22 
&23 



■1 + C 

4 + 2° 
1^2 



u tH -2C + C" 

j^TT 2 — |C + ^C 2 



-0.422... 
-0.461... 
-0.466... 
0.227... 



(4.22) 



where C = 0.577215... denotes Euler's constant. Combining Eqs. ( 4.17| ) and 
( |4.20f ) we find for M(t) an asymptotic expansion in inverse powers of logt, 



M(t) 



Tit 



log 8t 
irt 

+ 



m 



-fei 2 m + ivmi] 



log 2 8t 

+ ; — q— [&22"io - 27r6i 2 mi + 7r 2 m 2 ] 



log 3 8* 
irt 

log 4 8t 



[-6 32 m + 37r6 22 mi - 37r 2 6 12 m 2 + vr 3 m 3 ] 



(4.23) 



where we have set h 



'02 



1. We recall that whereas the b m £ are numerical 



coefficients, the m n defined by Eq. ( 4,18| ) are specific for the observable M. 
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Eq. (4.23) shows that the leading asymptotic behavior is ~ t/log8t for ob- 
servables that have mo ^ 0, and ~ i/log 2 8i for those that have mo = 
but mi ^ 0. In view of the definition of mo and the discussion at the end 
of Sec. |3.2| , observables with mo 7^ involve only black patterns and will be 
called, for short, "black" observables, whereas those with mo = will be 
called "black-and-white" observables. 



4.4 Examples 

It is now easy to derive results for many examples of interest by applying 
the formulas of Sec.fO. 



Example 1.1. We take for M the total number S of sites in the support. 
The first column of coefficients ha m Table I shows that the only nonzero 
term in the sum ( |4.18| ) is due to the set A = {0}. Since this set has ha = — 1 
and gA = 0, the only nonzero coefficient produced by Eq. ( [4.18 ) is mo = 1. 
From Eq. (4.23) we then have 



S(t)~ 



TTt 



1 



n-2 



+ 



>22 



>12 



log8U log8t log 2 8i log d 8t 



+ 



(4.24) 



The numerical values of the coefficients of the first two subleading terms are 
given m Eq. (Egg ) and agree with those of Torney pl| . 



Example 1.2. Next we take for M the total boundary length E between the 
white and black areas. In the second column of coefficients ha i n Table I 
shows that only two sets A enter, with the pair (ha, 9 a) equal to (4,0) and 
(-4, i). Eq. ( J4.18D then leads to m = and m n = 2 2 ~ n for n = 1,2,..., 
after which Eq. (4.23) gives 

2vr 2 t r 



E(t) 

with the coefficients 

£1 

So 



1 



£l 



+ 



£2 



log 2 8t L log 8t log 2 8t 



+ 



(4.25) 



Itt-2612 = 2.41636... 

I- 2 



}Trb 12 + 36 2 2 = 3.06116. 



(4.26) 



Example 1.3. Now take for M the total number I of islands. Using the 
third column of coefficients ha in Table I as input in Eq. ( 4,18| ) we find that 
mo = and 



m r , 



2 1 -" - [(7r + 2)/27r]" 
for n = 1, 2, . . . . Substituting as in the previous examples we find 



(4.27) 



j( t )~ ivr(vr-2) 



log 2 8t 



1 + 



log 8t 



+ 



+ ■ 



(4.28) 
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Analytic expressions for the coefficients ii are easily found with the aid of 
the preceding formulas but are of little interest. The numerical values of the 
first two of them are 



fci = - 2.087... 

t a = -21.304... (4.29) 

In the leading asymptotic behavior of both E(t) and I(t) an extra factor 
l/log8i appears compared to that of S(t) as a consequence of mo being 



zero, i.e., of E(t) and I(t) being black-and-white observables. In Sec. 4.5 
we shall come back to lit) and also make a comparison with the numerical 
simulations by Coutinho et al. Q. 

Example 1.4- Let (3 be a specific type of island, where type indicates shape, 
size and orientation. Let I pit) be the observable that counts the total num- 
ber of islands of that type. According to the preceding discussion we must 
have that 

Ip{t) ~ ivr(vr - 2)/^^ as t -> oo (4.30) 

for some proportionality constant fg, even though we cannot expect the 
approach to this asymptotic behavior to be uniform in (3. By summing 
Eq. flOcD on all (3 and comparing to Eq. ( p8|) we conclude that the average 



number of islands of type (3 represents, as t — > oo, a fixed fraction fp of 
the average total number of islands. We have calculated the fraction f\ of 
islands that are single isolated sites and the fraction /2 of "dimer" islands, 
consisting of two neighboring sites, with the result 

7r(7T 3 - 7ir 2 + 147r - 4) 

(7T — 1J(7T Z — D7T + 4) 

/ a = 0.073557... (4.31) 

The main effort goes into the calculation of the necessary coefficients qa- 
Those needed for I\ are listed in Table I. We do not present the 22 coef- 
ficients needed for I2, nor the final analytic expression for f2- Since the 
dimers may have two orientations, the islands of sizes 1 and 2 represent a 
fraction /1 + 2/ 2 = 0.707193... of all islands. 



With these results in hand we return to the total boundary length 
considered in Example 1.2. On an infinite lattice one can write E(t) = 
E ex tit) +Ei nt (t), where E ext (t) is the external perimeter of the support and 
Eintit) the total perimeter of the islands enclosed by it. In Example 1.2. we 
determined only the average of their sum; determining E ex t{t) and Ei n t(t) 
separately is a much more difficult problem that we have not seen how to 
solve by the present method. A rigorous lower bound for Ei n t(t) is neverthe- 
less easily obtained by adding up the perimeters of the single-site and dimer 
islands, and taking into account that all other islands have a perimeter at 
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least equal to 8. This gives Ei nt (t) > 0.4965. ..E(t). By arguments different 
from ours Lawler |^2| shows that in fact for t — ► oo the external perimeter 
increases only as E ex t(t) ~ i 2 / 3 , so that Ei n t(t) ~ E(t). 

Example 1.5. Let finally M{i) be equal to the pattern number N a (t) ob- 
tained by summing the expression over all sites x. If neither A± nor A2 
is empty then the pattern a is composed of both white and black sites and 
has mo = 0. Hence the average total number of these patterns is obtained 
by setting mo = in Eq. (4.23) and increases with t as 

N a {t)^ mi -^- (4.32) 
log 8t 

where mi is a dependent. The black-and-white patterns considered here 



necessarily lie on the boundary of the support; comparison of Eqs. ( 4.32 ) and 
(|4.25|) shows that the number per unit of boundary length of such patterns 
tends to a fixed value as t — > 00. 

4.5 Long-time behavior of averages. Finite lattice 

For a finite lattice of N sites, after an initial increase with time identical 
to what happens on the infinite lattice, we must expect deviations from 
the infinite lattice behavior to appear on a characteristic iV dependent time 
scale t(N) that will tend to infinity when N does. For larger times all black 
observables will level off and tend to a constant times N, and all black-and- 
white observables will pass through a maximum value, then bend down and 
asymptotically approach zero. The expansion that we shall look for in the 
case of a finite lattice will therefore involve a determination of this time 
scale. 

We shall not attempt in this case a full asymptotic expansion as for the 
infinite lattice, but only determine the leading asymptotic behavior. Our 
starting point is again Eq. ( 4.16j ), in which, when z — > 1, using Eqs. (|2.5|) 
and (2.6), we may substitute 

G (z) - g A (z) = —± + - log cN - g A + ■■ ■ (4.33) 

N{1 — Z) 7T 

where as before g A stands for <m(1) and the dots denote higher order terms. 
As for the infinite lattice, our procedure will be to bring the summation on 
A outside the integration on z and shift the integration path around the 
poles of the integrand. Therefore we should now discuss these poles. All 
required knowledge about the behavior of the various quantities involved has 



been collected in Sec. 2u. First, it follows from Eq. ( 4.33 ) and the definition 



flOp of m that for z -> 1 



E TT3 - -Nm Q {l - z) (4.34) 



A+% L Go( z )-9a(z) 
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Hence the integrand of Eq. ( |4.16| ) has a simple pole at z = 1 with residue 
Nmo. Secondly, as is clear from the discussion in Sec. 2^, this integrand 
has special simple poles for z = za, where 



ZA 



1 + 



7T 



N log cN 



1 + 



"^9 A 
log cN 



1 



log 2 iV 



(4.35) 



These are the only poles at a distance of 0(l/N\ogN) from z = 1, all the 
other ones being at least at distances of 0(1/N). Therefore, on time scales 
that are at least of 0(N log N), the other poles will contribute vanishingly 
to the result. 

Carrying the integral out but retaining only the contribution of the poles 
at z = 1 and z = za leads to 



M(t) ~ Nm 

+ N VA exp 



TTt 



N log cN 



1 + 



^9 A 
log cN 



+ 0{ 



\og 2 N 



(4.36) 



The first term in this equation is not present for the black-and-white ob- 
servables, which have mo = 0. We shall discuss these observables first. Two 
different time scales are of interest. 



1. The main regime, in which an observable takes values of the same order 
as its maximum value. To focus on this regime we scale time as 



TTt 



N log cN 



(4.37) 



and take the limit N ^oo, t^ooatr>0 fixed. In this "r -limit" Eq. ( |4.36| ) 
leads directly to 



77/ \ KN 

M (t) ~ — mi r e 



log cN 



(4.38) 



with mi defined by Eq. (fDJ|). 



2. The long-time regime, in which the observable approaches its final value. 
It is now appropriate to scale time as 



Tit 



a 



Nlog z cN 

with a > 0. In this "a -limit" Eq. (p~36| ) gives 

M(t) ~ N(cN)- a VAe~ ngA(J 
A^% 



(4.39) 



(4.40) 



One may notice that the result ( [4.38 ) is recovered from Eq. ( 4.40f ) if one sets 
a = t/ log cN and expands the terms inside the sum on A to first order in 
1/log cN. 
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For black observables Eq. ( 4.36Q with the scaling of Eq. ( 4.371) leads di- 
rectly to 

M(t) ~ Nm (l - e~ r ) (4.41) 

This decay law depends on the type of the black pattern only through the 
prefactor rriQ. 

4.6 Examples 



We consider the same examples as in Sec. 4.4 but now on the finite lattice. 



Example 2.1. Let M = S. Substituting mo = 1 in Eq. ( 4.36 ) we see that on 



a finite lattice of N sites the average number of unvisited sites decays as 

N - S(t) ~ Ne~ T (4.42) 



a result first obtained by Weiss et al. |24| . It is also valid in the long-time 
regime, where it can be written as N(cN)~ a . In this regime, since a > 0, 
the unvisited sites constitute an infinitesimally small fraction of all lattice 
sites. 



Example 2.2. Let now M = E. Using in Eq. ( 4.36Q the appropriate coeffi- 
cients ha from Table I we obtain 

E(t) ~ T- T Te~ T 

v ' logciV 

E(t) ~ AN(cNy a (l - e-™ 12 ) (4.43) 

in the main and long-time regimes, respectively. The prefactor N{cN)~ a in 
the second one of these equations is the average number of unvisited sites 
found in the preceding example. If these sites were randomly distributed, 
then since they are infinitely dilute, each of them would have four visited 
neighbors and the result for E(t) would be only 4N(cN)~ a . Hence the factor 
1— e _7r<7 / 2 represents nontrivial correlations due to unvisited sites clustering 
together. 

Example 2.3. Let M = I. Using in Eq. fl4.36| ) the coefficients /ia from Table 
I appropriate to this case, we find 

v ; logciV 

7(t) ~ N(cN)- a (l - 2e-™ /2 + e -^+ 2 W 2 ) (4.44) 

for the main regime and the long-time regime, respectively. Again, the factor 
in parentheses in the last equation is due to correlations in the positions of 
the unvisited sites. 

This example requires the following remark. The expression of Eq. ( |3.3| ) 
for the number of islands I(t) is correct only for the infinite lattice. For 
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the finite lattice with periodic boundary conditions in both directions one 
should add —1 when the support closes onto itself around the torus in one 
of the directions and —2 when it does so in both directions. In the latter 
case the first term in Eq. (|3.3|) is absent and one obtains the expression that 



we used to derive Eqs. (4.28|) and (4.44). The extra term —2 is of course of 



no importance in the main regime, but it needs to be taken into account in 
the long-time regime to ensure that I(t) vanishes when t tends to infinity. 

If we neglect the difference between averages of ratios and ratios of av- 
erages, then we have from Eqs. ( 4.42| ), (4.43), and (4.44) that the quantities 



[N - S(t)]/I(t) ~ 1/(1 - 2e" mT/2 + e -^+ 2 W 2 ) 

E(t)/l(t) ~ 4(1 - e-^ CT/2 )/(l - 2e~ mT/2 + e -( w + 2 W 2 ) 

(4.45) 

represent the average area and the average perimeter, respectively, of an 
island. In the main time regime the expression for the average area simplifies 
to 

[N-m/m^^fr- 1 (4.46) 

It is now possible to make a comparison with the simulations by Coutinho 
et al. U] , carried out on finite square lattices of up to N = 1200 2 sites 
with periodic boundary conditions. These authors were interested in the 
"fragmentation" of the finite lattice into islands and the way the average 
number and size of the islands eventually tend to zero. The comparison 
leads to the following conclusions. 

(i) In an early time regime, which for a lattice of 600 2 sites corresponds 
to t less than ~ 0.5 x 10 6 , the finite lattice size still plays no role and our 
Eq. (gjp for I(t) agrees within error bars with the simulation data shown in 
Ref. [||] for N = 600 2 . These data are for t > 0.05 x 10 6 ; for the agreement to 
be of this quality, not only the leading order term but also the two correction 



terms in Eq. ( 4.28| ) have to be taken into account. 



(ii) In the early time regime and in the main regime, where I(t) passes 
through its maximum, our long time expansion Eq. (4.44) overestimates the 
numerical values of I(t) by up to 50%; in the long time regime, when I(t) 
starts to decay, the agreement with the simulation data becomes rapidly 
better and stays good all the way up to t ~ 19 x 10 6 , where with a large 
probability no unvisited sites are left. 

(Hi) Coutinho et al. in their simulation find the time dependent average 
island size on a lattice of N sites to be a function only of the scaling variable 
a. Our expression (4.45) confirms this result. The scaling function is not, 
however, the power law that it was thought to be in Ref. M , but the inverse 



of a sum of exponentials given in Eq. (4.45). According to our Eq. ( 4.46 ) a 
power law appears only in the main regime and has the exponent —1. In 
the long time regime Eq. (4.45) leads to an apparent exponent with larger 
absolute value. 
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(iv) In the preceding discussion we have compared the numerical data 
for the average island size [N — S(t)]/I(t) to the theoretical result (4.45) for 
[N — S(t)]/I(t) that we were able to calculate. We have not been able to 
calculate directly the average island size. Nor have we been able to calculate 
still another quantity determined in the simulation of Ref . Q , viz. the time 
dependent "diversity" of the island sizes, defined as the number of different 
sizes that occur at any given time. 



5 Fluctuations and correlations 

5.1 Relation to a first passage time problem 

The original determination by Dvoretzky and Erdos || of the average num- 
ber S(t) of lattice sites in the support of a random walk was followed only 
much later [p], 11] by a calculation of the root-mean-square deviation of 



this quantity from its average. Yet that calculation was important, because 
the result, exhibited in our Eq. ( |1.2j ), shows that in the limit t —* oo the 
probability distribution of S(t) becomes infinitely narrow, even though only 
logarithmically slowly with t. 

It is now natural to ask if the more general observables whose averages 
we studied in Sec.|I| also have infinitely narrow distributions for t — > oo. 
Without much extra effort it will be possible to calculate also the cross- 
correlations. We therefore consider two observables of the form (|3.6|), namely 



M(t) = ]T fi A M A (t), M'{t) = ]T n' B M B (t) (5.1) 
A B 

which have J2a /M = J2b Mb = Oj an d focus on 



AM(t)AM'(t) = M(t)M'(t) - M(t) M'{t) (5.2) 
In this section we confine our analysis to the infinite lattice. 



The first step will be to find an expression for the generating function 
Cmm' defined by 

oo 

Cmm> (z) = E z< M{t)M'(t) (5.3) 
t=o 



When working out the RHS of Eq. (|J) with the aid of Eqs. |y) and (gj) 
we encounter averages of products m x+J 4(i)m y+ £(t) where x and y are ar- 
bitrary lattice vectors. It is then convenient to write y = x + r and to define 
U(r) as the union of A and r + B. 

A and B are finite sets, small in most applications, and we may locate 
them near the origin without loss of generality. For r not too small, therefore, 
U(r) will be the union of disjoint sets A and r+B, and its number of elements 
will be \U\ = \A\ + \B\. 
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With this notation we get 



M{t)M'{t) = Y J ^' B Y. M U(r){t) (5.4) 

A,B r 



The observable Mjjr r \ that occurs here is exactly of the form (|3.5| ), with A 
replaced by U(r). 

The generating function Cmm'( z ) is now easily expressed in terms of the 
first passage probabilities defined in Sec. |3.2| . The derivation is analogous to 
the one of Eq. (|4,5|) and the result is 

Cmm'(z) = -t— 2*MAA»b£ 2 Fu(r)(u,z) 

Z A,B r U £(7(r) 

= H VA» B J2[ E ^(r)( U ^) 

- y EF A (*,z)-Y l FB0>,*j\ (5-5) 

aGA heB 

in which the asterisk indicates that the term with A = B = is excluded 
from the summation. The function i*j/( r ) ( u i z ) i n Eq. ( |5.5| ) should be deter- 
mined from the linear system of equations 

J2 F u{r) (u',z)G(u-u',z) = -}— (5.6) 
u'£C/(r) z 

for all u 6 f (r). Formally this equation is strictly analogous to Eq. (4.7), 
but the structure of the matrix involved is different. It depends on the 
parameter r, which represents the distance between the two components of 
the set U(r). 

5.2 Scaling limit and long time behavior 

We shall be able to solve -Ft/(r) from Eq. ( f>.(\ ) only in the scaling limit z — ► 1, 
r — ► oo, with £ 2 = 4r 2 (l — z) fixed. But since z — ► 1 is exactly the limit of 
interest, this solution suffices provided the sum on r in Eq. fl5.4|) is dominated 
by values r ~ (1 — z)^ 1 / 2 . In the scaling limit we need to consider only the 
case where A and r + B are disjoint. Using the shorthand notation of 
Eq. (p| we find from Eq. (O) 



( X ~ 7aa'(2))-?V + 7a,r+b'0)-7 r r+b' = 1 

a'GA b'eB 

53 7r+b,a'0s)-7v + ^ (1 -7bb'(^))^r+b' = 1 (5-7) 



a'GA b'eB 



valid for all a £ ^4 and b £ B, with J 7 now defined on the set £7(r). In the 
scaling limit Eqs. (|2.14| ) and ( f4.8|) imply the behavior 

7a,r+bW ^ 2TT- 1 K (2rVT^)/G (z) + 0(1 - z) (5.8) 
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whose important feature is that to leading order no dependence on a or b 
appears. Furthermore 

, , q(a — a') , , ofb — b') 

7aa'(*)- ^ 7- / , 7bb'(*)-^W^ (5-9) 

and we shall henceforth suppress the argument of the 7's. 

In the scaling limit Eqs. (|5.7|) can then be written, in matrix notation, 

as 



j[a,a] _ 7 (A) \jl a >P] 

Xj\PM jiM _ 7 (B) 



^ = j (5.10) 

where a = (3 = \B\, 7^) and 7^ are matrices as defined in Sec. 4^ 



j[^. m l i s th e £ x m matrix with all elements equal to 1, j is the vector with 
all elements equal to 1, and 

\(Z,z) = 2ir- 1 K (a)/G (z) (5.11) 

We can obtain from Eq. (|5.10| ) a formal expression, analogous to Eq. (PH) , 
for the sum of the components of T. In the appendix it is shown that this 
expression can be transformed to 

V- t 2 - 2A - 7a - 7B / K10 x 

2^ -^u = 1 — T2 : 5.12 



Upon substituting Eq. ( 5. 12j ) in Eq. fl5.5|) and replacing the sum on r that 



occurs there by an integral on 4 we find 

Cmm 



2(1 - zYGo(z) A B ^ Jo (5 13) 



in which the function 2) is given by 
2 - 2A - 1A - IB 



1 - A 2 - JA - IB + 1A1B 1 - 7A 1 - IB 



(5.14) 



It is useful to observe that the three quantities ja, 7_b, and A(£, z) are all of 
order Gq 1 (z). 



The steps that follow are again analogous to the procedure of Sec. |4.3| . We 
wish to expand I(£,z) in inverse powers of Gq{z) and substitute the result 
in Eq. ( 5,13| ). The fi' B that characterize the observable M' define coeffi- 



cients m' n analogous to the m n of Eq. ( |4.18| ) . In the expansion we encounter 
furthermore the coefficients 



OO 







d££KS(0 (5.15) 



of which we shall need the explicit values a\ = 1 and = \, as well as 

1 r 1 log? 

°3 = -9 / ^ 1 TT72 = 0-58597... (5.16) 

2 Jo 1—4+4 
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Finally, before working out this expansion it is useful to classify the ob- 
servables M according to their order. We shall say that M is of order k 
if 



uiq = mi 



m k-i = and m k ^ 



(5.17) 



Observables of order k = and k = 1 have occurred in the preceding 
sections, and physically interesting examples with k > 2 perhaps exist. Let 
now k and k' be the orders of M and M' , respectively. We anticipate - 
as will be confirmed by the calculation - that we have to expand Cmm'( z ) 
in Eq. ( pU3| ) to order 1 / G^ +k ' +4 {z) . In view of Eqs. ( fUTp and ( gig ) only 
those terms in the expansion will survive the summation on A and B that 
contain at least a factor 7^ and a factor 7^ . Since there is one factor 



1/Gq(z) outside the sum in Eq. ( 5.13 ), this leaves room for at most three 
factors A. We have therefore found it convenient to begin by expanding 
/ = I'X + I"X 2 + /"'A 3 + • • • and then to determine the first three coefficients 
of this series in terms of 7,4 and 75. Then the 1/Gq{z) expansion of Cmm'( z ) 
leads to 



Cmm'(z) - 



1 



(1 



z) 3 G k+k' + 2 (2) 



2aimkm' k / 

-Gq l (z)(ir~ 1 (k + k' + 2)a 2 m k m' kl - 2ai{m k m' k , +l + m k+ im' k ,)^ 

+G 2 {z)(8Tr- 2 (k + l)(k' + l)a 3 m k m' k , 

— 2TT~ 1 (k + k' + 3)a2(mkm' k / +1 + m^im'j./) 
+ 2ai{m k m' k , +2 + m k+1 m' k , +1 + m k+2 m' k ,) 



+ 



(5.18) 



The next step is to invert Eq. ( p.3[ ). After integrating on z with the aid of 
Eq. flPcp one finds 



M(t)M'(t) ~ 



7T 



fc+fc'+2^2 



log 



8/ 



a\m k m k i 



log 8t 
1 

log 2 8i 



+-^— (-(A; + fc' + 2)(2ai6i 3 + a 2 )m k m' k , + 7rai(mi.m' fe / +1 + 
loe 8r v 

fc + fc' + 2){k + k' + 3)(ai6 23 + 2a 2 fei 3 ) 
+ 4(* + l)(k' + l)a 3 )}m fe m / i ., 
- (A; + k' + 3)7r(o 2 + 2aibi^){m k m' kl+l + m k+ im' k ,) 
+ 7r 2 ai(m fc m^/ +2 + m k+1 m' k , +1 + m fc+2 m' t ,)j 



+ • 



(5.19) 
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From Eq. (4.23) we deduce that when M is of order k its average is given by 

M{t) ~ \og k+1 m \- mk + ioTst V " ( k + 1)6l2mfe + nmk+1 



+ i — Y7^(k( k + + 2 )b22m k - ir(k + 2)bnm k+ i + TT 2 m k+2 
log 8t v 



+ 



(5.20) 



Eq. (5.20), its counterpart for M'(t), and Eq. ( 5.19| ) now have to be combined 
in Eq. (|5.2j ). Using the explicit expressions for the coefficients a\, 02, &12, 
622 j an d 623 leads to the desired correlation function. For t — > 00 the two 



leading orders cancel in the subtraction in Eq. ( |5.2|) . The final result is 



1 AM(t)AAf(t) A^ 

(fc + l)(fc' + l) M{t)M'{t) log e VU l ' ; 

where the dots indicate terms of higher order in 1/ log 8t and 

A e = AH 3 + 00 - ^vr e = oo.^VBi.. (5.22) 

Eq. ( |5.21| ) contains as a special case the well-known result of Eq. ( |1.2| ) for 



the variance of S(t), originally due to Jain and Pruitt [1C], and rederived 



with the aid of a method more similar to ours by Torney flllj . 



5.3 Conclusions 



It is remarkable that the ratio on the RHS of Eq. ( [5.21 ) is universal: It is 



independent of the choice of the observables M and M' . But we shall now 
see that Eq. (]5.2l| ) has consequences that reach far beyond this simple fact. 

For k = 0, 1, 2, . . . we define for each observable M of order k its nor- 
malized deviation from average, tjm, by 

log8t AM(t) 
m{t) = (k + l)A M(t) (5 ' 23) 

These variables satisfy to leading order 



VM(t) = 0, VM(t)r)M'(t) = l for all M, M' (5.24) 

It follows that for any two M and M ' the difference tjm (t) — tjm' (t) has zero 
variance, and therefore the r]M(t) are all equal to a single random variable 
that we shall call rj(t). As a consequence we can relate the deviation from 
average AM(t) of any observable M(t) of order k to rj(t) by 

AM(t) ~ (k + M(t) r!(t) (5.25) 

log VU 

Upon writing down this equation for the special case M = S, using the 
explicit expression ( [4.24j ) for S, and eliminating rj(t), one finds 

is 

AM(t) ~ (jfe + 1) \ m k AS(t) (5.26) 
log'' 8t 
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This last equation embodies one of the main conclusions of this work: All the 
observables M fluctuate around their averages in strict proportionality with 
the fluctuation of total number of sites S(t) in the support. This conclusion 
applies, in particular, to the pattern numbers N a (t), the total perimeter 
length E{t) of the support, the total number I(t) of islands enclosed by it, 
and the total number Ipif) of islands of a specific type (3. We are not aware 
of any computer simulations that confirm Eq. ( 5.26| ), although they would 
be easy to carry out. 

There is a still different and instructive way to formulate this conclusion. 
Let M be an observable of order k and p a quantity that remains of 0(oo) 
when t — > oo. Eq. (5.20) can now be used to establish the asymptotic expan- 
sion in powers of l/log8t of p~ 1 M(pt). Upon comparison with Eq. ( |5.25f) 
and choosing log p = — An(\J) one finds that all fluctuating observables M(t) 
can be written to second order in the form 

M{t) ~ e A ^ u W(e- Ari{u) t) (5.27) 

In the mathematical literature on Brownian motion (for a review see Le 
Gall [jig|] ) the quantity —Atj/^ti) has appeared in the study of the asymp- 
totic behavior of the volume of the Wiener sausage (where it is commonly 
denoted by the symbol 7) and is known as the renormalized local time of 
self-intersections, a concept introduced by Varadhan in an appendix to an 



article by Symanzik [29]. 



6 Discussion 

In a preceding letter |J a more general approach was presented to the much 
more restricted problem of how to calculate the average I(t). The total 
number of islands was written as I(t) = C{t) — D(t), where the increments 
AC(t) = C{t) - C(t - 1) and AD(t) = D(t) - D(t - 1) are the numbers 
of islands created and destroyed, respectively, in the tth step (either AC(t) 
or AD(t) vanishes). Hence C(t) and D(t) are the total number of islands 
created and destroyed, respectively, up until time t. It was shown || that 
asymptotically for t — > 00 to leading order 

C {t) ^D {t) ^A^ t (6.1) 

where A = 0.1017... In the difference 7(t) = C(t) - D(t) the leading or- 
der ( |Q| ) cancels and the result fl4.28|) appears. It does not seem possible 
to express C(t) and D(t) as observables of the type M(t). The additional 
determination of C(t) and D(t) makes the calculation of Ref. || more in- 
volved. Nevertheless, the generating function found there (Eq. (5) of Ref. ||) 
for d 2 I(t)/dt 2 is equivalent to the one for I(t) implied by Eq. ( |4.16 ) of this 
work. 
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The success of the study presented here is due to the generating function 
method, whose potential is fully exploited. We also run into what may be 
the limitations of this method. There are several quantities that appear 
naturally but whose averages and variances we do not see a way to determine. 



These include the total external perimeter discussed in Sec. 4.4 and the area 
of the islands enclosed by the support. 

Many of the quantities discussed in the preceding chapters have close 
analogs in planar Brownian motion. We shall denote these analogs by the 
superscript B. The Brownian motion analog of the support of the lattice 
random walk is the set of points <Sjf (i) C that has been swept out in the 
time interval [0, t] by a disk of radius b performing Brownian motion with 
diffusion constant D. The set S B (t) is commonly called the Wiener sausage 
associated with the Brownian motion trajectory. The total area S^(t) of 
this set is analogous to the number of sites S(t) in the support of the lattice 
random walk. 

The diffusion constant D is defined so that the mean square displacement 
equals ADt. It can be scaled away, but we shall keep it here to facilitate 
comparisons between results from different sources. In the mathematical 
literature one customarily sets D = ^, whereas the long-time, large-distance 
limit of the random walk of this work yields D = j . 

The asymptotic behavior of S^(t) was first determined by Leontovitsh 



and Kolmogorov 16]. Berezhkovskii et al. give the complete asymptotic 
expansion 

b { > - log KDtb- 2 ^ log m KDtb~ 2 [ ' 

° m=0 ° 



where k = 4e 2C and the b m 2 are as defined by Eq. (4.21). The difference 



ASfr(t) was shown (T^] to be a random variable such that the distribution 
of t^ 1 (logt) 2 ASfr (t) converges for t — > oo to a limit distribution identical to 
the one of i~ x (log t) 2 AS{t). 

Two remarks can be made about the relations obtained by differentiating 
Eq. fl6.2|) with respect to b. First, let E^(t) be the total boundary length 
of the Wiener sausage S^(t). Assuming that dS^ is sufficiently regular we 
have 

E b B (t) = ^ (6.3) 
Upon averaging, using Eq. (|6.2|) and setting 602 = 1, we find from Eq. (]6] 



E ' ( '^S?^[ 1 + -] (6 ' 4) 

This expression has the same asymptotic time dependence as Eq. (|4.25| ) for 
E(t), but the coefficients do not coincide. Secondly, differentiate once more 
and consider the dimensionless number 

J b B (t) - (6.5) 
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Its average behaves in the large t limit as 



—s. , 8-n-Dtb- 2 , . 

J * (!) * (6 ' 6) 

It is easy to show that if the boundary dS^(\J) is sufficiently regular (having 
at least a tangent vector in each point), then the quantity ([T^) is equal to 
2ir(lf(t) — 2), with lf(t) the number of islands. However, these regularity 
conditions are not satisfied here since the boundary has cusp points. Never- 
theless, comparison of Eqs. and flPip shows that Jf(t) has the same 
asymptotic time dependence as I(t) (but with a different coefficient).. 

We now discuss the relation of our work to results that have appeared in 
the mathematical literature. Throughout the comparison it should be borne 
in mind that whereas those results are rigorous, the ones of this paper have 
been obtained by the usual methods of mathematical physics. 

Mountford (2(J was the first to study the connected components of the 
complement of the Brownian path of a point (i.e., the case b = 0). For all 
times t > the number of these components is infinite due to the presence 
of many small ones. However, one can ask for example what the number 
C £ (t) is of connected components with an area larger than a prescribed value 
ire 2 . Le Gall, strengthening the results due to Mountford, has shown that 
for almost all Brownian motion trajectories (with D = |) in a fixed time 
interval 

e 2 log 2 Dte- 2 
lim S CM) = 2 (6.7) 

e^o 4Dt w v ' 

where we have used dimensional analysis to restore the variables D and 
t. If one now assumes that C £ (t) is of the same order of magnitude as 
the total number if it) of islands in the Wiener sausage associated with 
the same Brownian motion trajectory executed by a disk of radius e, then 
Eq. ( |6.7| ) has an asymptotic time dependence that agrees with the result of 
our Eq. (|4.28|) . 



Werner |2l| considers the shape of the connected component containing 
a prescribed point not on the trajectory, and shows that for t — > oo the prob- 
ability distribution of this quantity tends to a well-defined limit law. This 
result is possible only because of the scale invariance of Brownian motion 
and it would be cumbersome to formulate its lattice counterpart, even in the 
long-time, large-distance limit. The statement of our Eq. ( |4.30D about the 
number of islands I pit) of given type /3, and the associated result about its 
fluctuation AIp(t) implied by the discussion at the end of Sec. 5, constitute 
the point of closest approach between this paper and Werner's. 

Finally, we summarize the new results of this paper. On the one hand, we 
give the explicit asymptotic behavior as t — > oo for several new observables 
associated with the support of the lattice random walk. These include the 
total perimeter length, the total number of islands, and the total number of 
single-site and dimer islands. On an infinite lattice the asymptotic behaviors 
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all consist of a leading term multiplied by a series in inverse powers of 
1 / log 8t. The perimeter and the number of islands are also considered on a 
finite lattice, where scaling laws are obtained in terms of the time t and the 
lattice size N, and a comparison with computer simulations by Coutinho et 
al. Q is possible. 

On the other hand, there is the important general result of Sec. pT3| , 
stating that the pattern numbers all fluctuate in strict proportionality with 
one another and with the total number S(t) of sites in the support. The 
fundamental fluctuating variable, called r) in this work, is the renormalized 
local time of self-intersections. This result strongly contributes to shape our 
picture of the support of the two-dimensional random walk: At any fixed 
time t, a large class of detailed properties is determined by the value of the 
single random variable rj. 

One of the open questions that can now be formulated is connected ex- 
actly with this random variable, which appears in our work as the time 
dependent variable rj(t). Existing results seem to concern exclusively its 
stationary distribution V{rf), which prevails in the limit t — > oo. Our in- 
vestigations point towards the interest of also studying the autocorrelation 
function rj(t)r](t') and possibly other time dependent properties. A second 
question that naturally comes up is: How does the picture in higher di- 
mensions differ from the one found here in d = 2? It should certainly be 
expected that in sufficiently high dimension the pattern numbers become 
independent random variables. The mechanism by which this independence 
comes about seems worthy of further elucidation. We shall leave these and 
other questions for future work. 
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A Appendix 

In this appendix we prove the matrix algebra results we use to derive aver- 
ages and correlations. 



The following result is used in Sec. [1^. Let A denote an invertible matrix 
of dimension I x £, with £ > 2. In what follows j[ m ' n l stands for the m x n 
matrix whose elements are jj™'^ = 1. Define furthermore 

9a 1 - J2( A '% = THJ^U- 1 ) (A.l) 
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We achieve our aim of finding a simplified form for 



(A.2) 



by rewriting Eq. 



as 



r A = Tv[j^U- 1 (i + jMa- 1 )- 1 



(A.3) 



where 1 is the unit matrix. 

An intermediate step of the demonstration consists in noting that 



j[k,l] CJ [£,rn] D = g-\j[*M D ( A .4) 

where C and D are square matrices of dimensions I x t and m x m, respec- 



tively. In Eq. ( A.3 ) we now expand the argument of the trace in powers of 
jMA' 1 , Iteration of Eq. flOl) and use of Eq. lead to 



Tr 



whence 



(j^U- 1 )" 



9a 



1+9A 



(A.5) 



(A.6) 



When I = 1 formula (A.6) remains valid for all A 6 M (and in particular for 
A = 0) if for that case one supplements Eq. ( |A.1| ) with qa = A. 

A generalized version of the above result is needed in Sec.||. It involves 
two invertible matrices A and B of dimensions i x i and mxm, respectively, 
with £, m > 2. We now wish to find a simplified expression for 



r^(A) = £ 

Upon using the decomposition 
/ + A AJ^ ,m l \ _ 



1 + jl^A- 1 xj^b- 1 



(A.7) 



A 
B 



(A. 



and writing the double sum in Eq. (A.7) as a trace with the aid of the matrix 
j[i +m ,e +m ] one findg 



where 



r AB {\) =Tr[M(l)(l + M(A))- 1 



J^A- 1 XJI^B- 1 



(A.9) 



(A.10) 
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We now expand the argument of the trace in Eq. (|A.9| ) in powers of 
M(A). Hence we are left with the calculation of Tr(M(l)M(A) n ), which is 
carried out by finding recursively in n relations between the blocks of 

M(l)M(X) n = U n =( X z ^ ^ (A.ll) 

Knowing that U n +\ = U n M(X), one has 

X n+1 =X n J^A~ 1 + XY n J^A~ l 
Y n+l =XX n J^B~ 1 + YnJ^'^B- 1 (A. 12) 



Hence, again with the use of Eq. ( A. 4 ), one proves recursively that X n and 
Y n are of the form 

x n =6,^-1 

Y n =r ?n j[ £ H jB -i (A . 13) 
with initial values Co = % = 1 an d the recursion relations 

Vn+i = A^ 1 ^ + g^rjn (A. 14) 

Upon setting Z n = £ n j[ m '^A -1 and T n = r n j[ m ' m l B" 1 one finds in a similar 
way 

Tn+l = gB l T n + A#^ Cn 

Cn+l = ^g B l T n + g^Cn (A.15) 

with initial values tq = Qq = 1. Hence, using that TrX n = g^ 1 ^ and 
TYT n = g]fT n , we have from Eqs. (|A.9| ) and (|A.10|) 

oo oo 
n=0 n=0 (A. 16) 

By summing each of the equations ( A.14| )-( [A,15| ) on all n with weight 
(— l) n one obtains two sets of linear equations from which X)n^=o( — l) n £n 
and X^o( — l)™ 7 "™ ma y be solved. Substitution in Eq. ( [A. 16 ) yields the final 
result 

r M X) = + (A. 17) 

i-X z +g A + gB+ gAgB 

By redoing the calculation with I or m = 1 one finds that this expression 
for remains valid for those special cases. 



30 



References 

[1] G. Polya, Math. Ann. 84 (1921) 149. 

[2] G. H. Weiss, Aspects and Applications of the Random Walk (North- 
Holland, Amsterdam 1994). 

[3] B. D. Hughes, Random Walks and Random Environments, Vol. 1: Ran- 
dom Walks (Clarendon Press, Oxford 1996). 

[4] K. R. Coutinho, M. D. Coutinho-Filho, M. A. F. Gomes, and A. M. 
Nemirovsky, Phys. Rev. Lett. 72 (1994) 3745. 

[5] S. Caser and H. J. Hilhorst, Phys. Rev. Lett, xx (1996) xxxx. 

[6] A. Dvoretzky and E. Erdos, 2nd Berkeley Sympos. Math. Stat, and 
Prob. (University of California Press, Berkeley 1951) p. 33. 

[7] G. H. Vineyard, J. Math. Phys. 4 (1963) 1191. 

[8] E. W. Montroll, Proc. Symp. Appl. Math. Am. Math. Soc. 16 (1964) 
193. 

[9] E. W. Montroll and G. H. Weiss, J. Math. Phys. 6 (1965) 167. 

[10] N. C. Jain and W. E. Pruitt, Z. Wahrscheinlichkeitstheorie verw. Ge- 
biete 16 (1970) 279. 

[11] D. C. Torney, J. Stat. Phys. 44 (1986) 49. 

[12] J. Rudnick and G. Gaspari, J. Phys. A 19 (1986) L 191. 

[13] S. J. Sciutto, J. Phys. A 28 (1995) 3667. 

[14] G. H. Weiss and S. Havlin, Phil. Mag. B 56 (1987) 941. 

[15] M. Kac and J. M. Luttinger, J. Math. Phys. 15 (1974) 183. 

[16] M. A. Leontovitsh and A. N. Kolmogorov, Phys. Z. Sowjet. 4 (1933) 1. 

[17] J.-F. Le Gall, Ann. Probab. 16 (1988) 991. 

[18] J.-F. Le Gall, in: Ecole d'ete de probabilites de Saint-Flour, I. Lecture 
Notes in Mathematics 1527, Ed. P.L. Hennequin (Springer, Berlin 
1992). 

[19] A. M. Berezhkovskii, Yu. A. Makhnovskii, and R. A. Suris, J. Stat. 
Phys. 57 (1989) 333. 

[20] T. S. Mountford, Stochastics 28 (1989) 177. 

[21] W. Werner, Probab. Theory Relat. Fields 98 (1994) 307. 

[22] G. Lawler, Elect. Comm. in Probab. 1 (1996) 29. 

31 



[23] W. Th. F. den Hollander and P. W. Kasteleyn, Physica A 112 (1982) 
523. 

[24] G. H. Weiss, S. Havlin, and A. Bunde, J. Stat. Phys. 40 (1985) 191. 

[25] M. J. A. M. Brummelhuis and H. J. Hilhorst, Physica A 176 (1991) 
387; 185 (1992) 35. 

[26] G. Zumofen and A. Blumen, J. Chem. Phys. 76 (1982) 3713. 

[27] F. Spitzer, Principles of Random Walk (Van Nostrand, Princeton N.J. 
1964). 

[28] F. S. Henyey and V. Seshadri, J. Chem. Phys. 76 (1982) 5530. 

[29] K. Symanzik in: Euclidean Quantum Field Theory, R. Jost ed. (Aca- 
demic Press, New York 1969). 



32 



FIGURE CAPTIONS 



Figure 1. A random walker makes t random steps between the centers of 
neighboring squares. The squares visited are colored black and constitute 
the support of the walk. In the example of this figure the support encloses 
four islands of unvisited sites. The boundary of the support is represented 
by a heavy line. 

Figure 2. Let the boundary of the support be oriented such that as one 
proceeds along it the white sites are on the left and the black ones on the 
right. Then in (a) the boundary turns through an angle +^vr and in (b) 
through an angle — ^ir. In (c) the black squares are (by convention) con- 
sidered to separate the white squares from one another and the part of the 
boundary shown makes two turns through 
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A 


S E I h 





{0} 

{<W 

{0,ei,e 2 ,ei + e 2 } 
{0,ei,e 2 } 
{0,ei,2ei} 
{0,ei,-ei,e 2 } 
{ei,e 2 , -ei, -e 2 } 


1 

-14 11 

-4 -2 -4 
1 

4 
2 
-4 
1 




1/2 

(tt + 2)/2tt 
7r/2(7r- 1) 
tt/4 

7r(7r-6)/2(7r 2 -6tt + 4) 
1 



Table I. The coefficients \ia for the four observables S, E, I and h defined in the 
text; entries not shown are zero. Symmetry under rotations over 7r/2 has been 
exploited to reduce the length of the table. The coefficients in each column add up 
to zero. The parameter gA is given for all sets A that occur; it is undefined for the 
empty set 0. A itself is defined only up to a translation. 
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